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, Abstract. We study the efficiency of quantum tomographic reconstruction where 

the system under investigation (quantum target) is indirectly monitored by looking 

' ^ ' at the state of a quantum probe that has been scattered off the target. In particular 

O , we focus on the state tomography of a qubit through a onc-dimcnsional scattering 

■j of a probe qubit, with a Heisenberg-type interaction. Via direct evaluation of the 

^H associated quantum Cramer-Rao bounds, we compare the accuracy efficiency that 

•3 one can get by adopting entanglement-assisted strategies with that achievable when 

^H entanglement resources are not available. Even though sub-shot noise accuracy levels 

'~^ are not attainable, we show that quantum correlations play a significant role in the 

p,i estimation. A comparison with the accuracy levels obtainable by direct estimation 

^ (not through a probe) of the quantum target is also performed. 
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1. Introduction 

The possibility of reconstructing the quantum state of a system via measurements 
(quantum state tomography, QST) is a central problem in quantum information 
theory pQ, which poses a series of fundamental questions related to the fact that 
the state itself is not directly observable and that each given quantum measurement 
typically reveals only partial information on the observed system. In recent years, a 
great deal of work has been devoted to this issue and many important features have 
been recognized, including the fact that having at disposal several copies (say M) of the 
initial state, collective measurements are more informative than individual ones, e.g., 
see Refs. [H El |3] and references therein. In abstract terms, QST ultimately reduces to 
the ability of estimating the set of continuos parameters which define the expansion of 
an unknown state with respect to a reference basis of operators (say the set of Pauli 
matrices for a two-level system, qubit). As such, its ultimate accuracy limits can be 
evaluated by exploiting some general results of quantum estimation theory [H El El HO] 
(more precisely, of a part of the theory which directly deals with the estimation of 
continuos parameters). 

Up to date, most of the works focused on scenarios where the system under 
investigation can be directly accessed, i.e., posing no constraint whatsoever to the 
physical operations one may perform on it. In this context, for instance, the ultimately 
accuracy limit obtainable in the tomographic reconstruction of a qubit initialized in an 
arbitrary (possibly) mixed state have been set in Ref. jS] , by computing the associated 
quantum Cramer- Rao (CR) bound [H [5], in terms of the quantum Fisher information 
(QFI) matrix of the problem (see below for the precise definitions). In this paper, 
instead, we address the problem from a slightly different perspective, which captures 
an important aspect of many realistic experimental situations. More precisely, along 
the line set in Ref. [H] , we consider the case in which the system of interest (from now 
on the target X), can only be addressed indirectly via measurements performed on a 
probe which has interacted with it. In our model the latter is described as a quantum 
system A characterized both by external (e.g., momentum/position) and internal (e.g., 
spin) degrees of freedom, which the experimentalists are allowed to prepare in any initial 
configuration (also the target possesses external and internal degrees of freedoms but, 
for the sake of simplicity, only these last are supposed to be unknown, the external 
degrees of freedom being assigned by fixing the position of the target system). The 
tomographic reconstruction then proceeds by letting A and X interact via a scattering 
process and by measuring the final state of the former (or at least a part of it which has 
been scattered along some preferred direction). Indeed, the whole setting is devised in 
order to mimic the basic features of a standard (Rutherford- like) scattering experiment, 
where one tries to reconstruct the properties of a target system by firing probe particles 
on it and by looking at the way they emerge from the process. For the sake of simplicity 
we will limit the analysis to the case of ID scattering processes, and describe the internal 
degrees of freedom of A and X as two-level (spin) systems (similar models have been 



Entanglement- assisted tomography of a quantum target 




Mp'^a) 



Figure 1. Scheme of a programmable quantum channel. Here the CPTL map ^^ is 
defined as an interaction with an external system X through a given (fixed) coupling 
represented by the operator Sxa- In our problem the channel to be reconstructed 
describes the ID scattering of a probe A off a target qubit X, with Heisenberg-type 
interaction. This process induces a modification of the probe initial state p^^ according 
to the target initial state px{v) identified by the coordinates v. 



recently analyzed to study entanglement generation [T2| fT3]). 

It is worth noticing that the problem we are considering admits also an 
interpretation in the context of quantum channel estimation theory (see Ref. [lOJ and 
references therein). This is the theory, sometimes identified with the name of quantum 
metrology, which studies the efficiency of those schemes designed to recover information 
not on a quantum system, but on a quantum channel (quantum process tomography, 
QPT). We remind that in quantum mechanics quantum channels represent the most 
general physical transformations and are fully described by assigning completely-positive 
trace-preserving linear (CPTL) mappings [5], which act on the density matrices of the 
system (in our case the probe A). In quantum metrology the mapping ^^ is assumed to 
belong to a family of transformations identified by a set of parameters v, whose values are 
unknown and which we wish to recover by preparing the system in some fiduciary initial 
state p™ and by measuring the corresponding output state transformed by the channel. 
In our case the quantum channel to be estimated is the one that induces a modification 
on the probe A via its interaction with the target X, while the v's correspond to the 
parameters that define the (unknown) state px (see Fig. [I]). In the jargon introduced 
in Ref. [TB], this transformation belongs to the special class of programmable channel^ 
A well-known fact is that, in general, if M is the number of tests we perform in order 
to recover the actual values of v (each test consisting of applying the same channel $„ 
to a new copy of probe A) the statistical scaling of the associated uncertainty can be 

I Explicitly, these channels can indeed be parametrized by assigning a fixed interaction with an external 
unknown system. In a seminal paper [^ Nielsen and Chuang proved that the family of all unitaries 
acting on n qubits is not programmable, since an n-qubit register can encode at most 2" distinct 
quantum operations. Notwithstanding the non-universality of programmable channels, the authors 
also pointed out the possibility to program the family of unitary operations probabilistically. In this 
context, an interesting model was proposed for the case of single qubits in [7^ and implemented for 
photonic qubits in [8]. 
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reduced from the "standard quantum limit" (SQL) (or "shot noise" in quantum optics) 
1/vM scahng, to the so-called "Heisenberg bound" 1/M scaling, by the introduction of 
suitable quantum correlations between initializations of the various copies of A [101 [E] . 
However, this is not the case if, as in our case, the channels under investigation are 
programmable [TB]. As a consequence in our model no sub-shot noise scaling in M of 
the accuracy should be expected. For this reason we will limit the analysis to those 
configurations in which the M tests on X are performed by preparing M copies of A in 
the same initial state. The QFI matrix approach [H EJ [9] will then be used to evaluate 
the associated accuracy, optimizing it with respect to the initial preparation of A (with 
respect to both its internal and external degrees of freedom) and comparing it with the 
results obtained in Ref. [3j for the case of direct estimation. Most importantly, we will 
also study an entanglement-assisted (EA) strategy, where each copy of A is initialized 
into an entangled state of A with an external ancilla system B, which does not interact 
directly with X (see Fig. |2]), showing a clear improvement in the performance of the 
estimation process with respect to the non-EA (NEA) strategy in line with the findings 
of Refs. [ia[T6]. 

The paper is organized as follows. In Sec. [2] we introduce our scattering model and 
briefly review some basic aspects of quantum estimation theory. Then, in Sees. |3] and |4] 
we will discuss two different strategies for the tomographic reconstruction of the target 
qubit, with and without the help of quantum correlations in the probe preparation. In 
particular, in Sec. [3] we first derive the exact expression for the QFI matrix in the case 
of an EA configuration, in which the probe is initialized in a maximally entangled state 
with the ancilla system, and discuss some applications of the result in the evaluation 
of some functionals of the state of the target (specifically its purity and its azimuthal 
angle). In Sec. Ill instead, we compare the EA and NEA cases by focusing on a special 
configuration in which the target state is characterized by a single unknown parameter. 
The paper ends with Sec. |5]by summarizing our results in the light of future perspectives. 
The more technical aspects of the derivation are presented in a couple of appendices. 

2. The model 

In this section we introduce the model and set up the notation. 

2.1. ID scattering of the probe A 

Suppose that a target qubit X is fixed at a given position x = on a line, and that its 
unknown state is described by the density matrix 

Px{v) = ^-^, 0<|«|<1, (1) 

with cr = (cr^', cr^, cr^) being the Pauli operators, and v = {vx,Vy,Vz) the 3D Bloch vector, 
which represents the set of parameters we wish to recover via QST. In Ref. |TTj it was 
shown that v [and hence px{v)] can be obtained from the transmission and reflection 
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probabilities of a probe qubit A scattered off X through a point-hke interaction, which 
couples the internal degrees of freedom of the two qubits via the following Heisenberg- 
type Hamiltonian 

2 

H=!^ + giax-CTA)SixA). (2) 

Here, m and pa are the mass and the momentum operators of A, g is a positive coupling 
constant, and ctj are the Pauli operators acting on the qubit J (= X,A). Specifically, 
in Ref. [11] A was assumed to be initialized into a known input state \k){k\ ® p™ and 
injected from the left of the line with momentum hk > 0. This kind of systems can 
be considered as models for a magnetic impurity spin embedded along a ID wire [T7], 
such as a semiconductor quantum wire [18] or a single- wall carbon nanotube [19]. For 
instance, an electron is sent through the ID wire as a probe particle, and its spin state 
is resolved after the scattering by spin-sensitive filters [20]. Alternatively, an electron 
populating the lowest subband of the ID wire can undergo scattering from two double 
quantum dots [2T] to which it is electrostatically coupled. The state tomography of 
Pxiv) then proceeded by solving the associated scattering problem and looking at the 
state of A which emerges either on the left (transmitted component) or on the right 
(reflected component), or both. Such states admit a simple expression in terms of the 
scattering matrix Sxa of the process defined by the unitary operator 



SxA = dk \k){k\ ^S*xA+ dk \-k){k\ ® S;,A, (3) 

with 

C*t C* _|_ C'^t or _ 1 c<*t or , c'^t ot _ n f /\\ 

'-'XA'-'XA "T '-'XA'-'XA — -"-J '-'XA'-'XA "r '-'XA'-'XA — ^- \^) 

Here the |/c)'s represent the momentum eigenstates of the probe A {hk being the 
associated eigenvalues), while the 4x4 matrices Sxa ^^^ Sxa define the spin-dependent 
scattering amplitudes associated with transmission and reflection events, respectively. 
They are functions of the probe wave number k and are given by 

S*/a = Oit,ri^) + /3t,r(fi)(crx " cta), (5) 

with fl being the dimensionless parameter 

mg 

n\k 

and 



^=^rrri. (6) 



,„. _ 1 - 2zfi _ -in 

"*^^) - (l-3.fi)(l + .fi)' ^*^^^ - (l-3.fi)(l + zfi)' ^^^ 

"-'"'= (i-3.n)(i + .n) - "'(") = "'(")■ («) 

For an explicit derivation of all these equations we refer the reader to [11] . 

Consider now an experimental setting in which an observer tries to reconstruct 
px{v) by merging incoherentl}l§] the data associated with the transmission events (to 

§ By incoherent merging of the transmitted and reflected data, we mean that no joint measurements 
are aflowed on the transmitted and reflected signals (a scenario which is reahstic if the rhs and Ihs 
detectors are located sufficiently apart from each other). 
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the right of the target) and the reflection events (to the left). The final state of A can 
be expressed as the tensor product of effective qutrit density operators 

+ Ttx{S^Upx{v) ® p'X]S^J^} ® |e*)A(e*|, (9) 

where Trx{---} is the partial trace over X, and \e*)A (|e^)yi) is a vector orthogonal 
to all the internal spin states of A, which represents the vacuum state associated with 
no particle reaching the rhs (Ihs) detector. The first term represents the contribution 
associated with a transmitted A reaching the rhs detector (vacuum on the Ihs), while 
the second represents the opposite one (i.e., spin on the Ihs and vacuum on the rhs). The 
detectors are assumed to have 100% efficiency and no particle creation or destruction 
by the target is possible, so that the number of particles is conserved: every incident 
particle emerges either from the left or from the right of the target. 

On the other hand, if the observer collects only transmitted particles, emerging 
from the rhs of the target, he will either see nothing {A being reflected by X) or see A 
emerging with the same momentum hk it had when entering the line but with a modified 
spin state due to the interaction with X. Such configuration is described by the density 
operator, obtained from p*^^{v) by tracing out the reflected case, namely, 

= TrxiS'xAPxiv) ® p'X]S'^^} + TT{S^^^[pxiv) ® f/X]S;c\}\e')A{e% (10) 

Notice that here 

TT{S^xA[Pxiv) ® p'X]Sx\} = 1 - ^^{S'xAlPxiv) p'X]S'^^} (11) 

is the reflection probability (i.e., the probability of no detection of A on the rhs of the 
line) . 

An analogous expression holds for the alternative experimental setting in which the 
observer only collects information of the signals emerging from the Ihs of the line. In 



this case Eq. (10) is replaced by 

P:,{v)=Tt,p'+''{v) 

= ^r{S'xA[Pxiv) ® p'X]S'^A\elA{e^\ + Tr x{S^^A[Pxiv) ® p'X]Sx\}- (12) 

The output density matrices p^ ^ y")-, associated with the three different 
experimental settings, are functions of the state px{v) of X. Therefore, one can 
acquire information on the latter by performing QST on the former. Furthermore, 
since p^ ^ (v) also depends on the input state p™ and on the input momentum hk 
of the probe A, one can try to optimize the resulting accuracy with respect to these 
parameters. 

2.2. Entanglement-assisted scheme 

An interesting variation of the previous schemes is obtained by considering the case 
in which A is prepared into a joint (possibly entangled) state p™^ with an ancilla 
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system B, that is not directly interacting with X (see Fig. |2]). Such EA configurations 
proved to be successful in boosting the efficiency of several quantum estimation [1^ and 
discrimination schemes |23j. 

Assuming that B sits at rest in the laboratory of the observer, the resulting states of 
AB emerging from the AX interaction can again be expressed in terms of the scattering 
matrix Sxa given before. Specifically Eqs. ^, (10) and (12) become 



+ Ttx{{S^xa®W[Px{v) ® p'XbKS'Ja^^} ^\e')A{e% (13) 

PAsiv) = Tir Pa's 

= Ttx{{SI,a ® Ib)[px{v) ® P^2b]{S'^a ® Ib)} 

+ TixaHS^xa ® W[Px{v) ® p%]{S'x\ ® Ib)} ® \e')A{e% (14) 

PABi^) = Tr,p!+/ 

= |eOA(e1 ® T^txaHS'xa ® Ib)[px(^) ® p!4"b]('^?a ® ^b)} 

+ Tr^{(^^^®l5)[px(«)®p!?B](^x^^®M}, (15) 

where Ib stands for the identity operator on B. 

2.3. The Cramer-Rao hound 

In the following sections we will compare the accuracy one can get by reconstructing 
Px{v) through the EA configurations via measurements on the output states pj^j^ i^) 



defined by the Eqs. (13), (14) and (15), with the corresponding accuracy one achieves 
with the NEA configurations associated with the output states p^ '' (v) of Eqs. (9), 



( 10 ) and ( 12 ) . Differently from [1 Ij but in line with the approach of [Sj |9] , such accuracies 
will be evaluated by computing the corresponding quantum CR bounds [U El Ej. 

We remind that the quantum CR theorem establishes a fundamental lower bound 
on the uncertainty of any estimation strategy devised to recover the three components 
of V. Speciffcally, assume that one has M copies of the state p{v) which encodes such 
parameters. (In our case, for the NEA setting p{v) is given by the states p^ *" (v), 
depending on whether the observer collects only transmitted data, reffected data, or 
both. Similarly for the EA setting, where p{v) is identiffed with Pj^^ ('^))- ^ generic 
estimation strategy consists in assigning a (possibly joint) POVM measurement on the 
M copies of the state, and a classical data processing scheme that starting from the 
measurement outcome produces an estimation v*^^^ = (t>^***,t>|''*,t>3'^*) of the parameters 
V. The uncertainty of the estimation can then be evaluated by means of the covariance 
matrix Cov[v]jk = {v'j^^ — Vj){vl^^ — Vk) {j,k G {x,y,z}) obtained by averaging the 
distances between the real value of the parameter v and their estimations. In this context 
the CR bound implies that, independently from the adopted POVM and classical data 
processing scheme, such matrix must verify the inequality 

CovM > IjH{v)-\ (16) 
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where H{v) is the QFI matrix of the encoding state [U |9], i.e., 

^ y^ {d,Pn){d,Pn) + 2Re V ^-^^^^^{^|J^\^,^P^^){^,^P^\^Pn). (17) 

Here, p{v) = J2n Pnli^n) {i'nl is the spectral decomposition of the encoding state p{v), 
while dj stands for the partial derivative with respect to the jth component Vj of the 



vector V. The jth diagonal element of the inequality ( 16 ) provides the CR bound for the 
variance associated with the accuracy in the estimation of the parameter Vj, for fixed 
values of the others, i.e.jj]] 

1 



Cov[t;],, = Var^.] = {vf - v.f > -^{H-% (18) 



(notice the 1/a/M SQL scaling of y/Var[fj]). It is worth observing that even though 



the bound ( 16 ) is not always attainable, the bound ( 18 ) is known to be asymptotically 
achievable for a sufficiently large M. More generally, Eq. (16) permits also to derive 
an accuracy bound on the variance of the estimation of any given function / = f{v) of 
the parameters t^ [H [9] . Specifically by re-parameterizing the problem with a new set of 
independent parameters v = {vi{v),V2{v),V3{v)), which include the quantity of interest 
as (say) the first element vi{v) = f{v), one gets, 

Var[/] > j^iH-')u, (19) 

where H = BHB^ is the QFI matrix of the parameters v obtained from H via the 
similarity transformation induced by the Jacobian matrix Bjk = dvk/dvj. As in the 



case of Eq. (18), for fixed values of V2 and v^ the inequality (19) establishes a bound on 
the accuracy reachable in the estimation of the function / (the bound being achievable 
for a sufficiently large M under the assumption that V2 and v^ are known a priori). 

As an application of the above construction, and for future reference, it is instructive 
to report the quantum CR bounds associated with a direct estimation of pxiv) which 
has been first computed in Ref. [3]. In this scenario the observer is assumed to have 
complete (not probe-mediated) access to the target state X, so that the encoding state 
p{v) introduced above coincides with the density matrix px{v). In this case the QFI 
matrix possesses a simple form in polar coordinates, 

{vx,Vy,Vz) = {r sin 9 COS (f),r sin 9 sin (l),r COS 9), (20) 



where it is diagonal. Indeed using Eq. (17) one finds. 



Cov[(r,e,0)]>^/f^-(r,e,0)-\ (21) 

where i3"**''"(r, 9, 0) is the QFI matrix given by 

/f dir(^^ Q^ ^) _ diag(c^'^ cf\ cf "■ sin^ 9), (22) 

II Indeed, one can easily verify that the quantity [{H~^)jj]~^ coincides with the QFI function 
associated with the estimation of the parameter Vj on the one- parameter family of states p{vj) obtained 
from p{v) when assigning fixed values to the other components of v. 
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Figure 2. (Color online) Strategies for the estimation of the target parameters. In the 
NEA strategy the probe A and the ancilla system B are in a separable state, and, as will 
be clarified in Sec. |4j this analysis can be performed by completely neglecting B. On 
the other hand, in the EA strategy quantum correlations are introduced (rectangular 
box). In both schemes one makes use of optimal POVMs for the estimation of the 
target parameters. Notice that only qubit A undergoes a direct interaction (Sxa) 
with the target X. 



dir 



^dir 



(23) 



Notice that the matrix i?'^"' does not depend upon the azimuthal angle while it is a 
function of the radial coordinate r and on the polar angle 6 (the latter however is just 
a geometric artifact introduced by the polar coordinates: the north and south poles are 
indeed insensitive to rotations along the 2;-axis). 



3. Entanglement-assisted strategy 

In this section we will consider the case in which the probe A is prepared in an entangled 
state with the ancilla B, that is the EA strategy introduced in the previous section and 
represented in Fig. |2j More precisely, we will assume that subsystem AB before the 
scattering is in a maximally entangled state, and compute the analytic expression for 
the associated QFI matrix i? as a function of the set of parameters defining the initial 
state of X and of the incident momentum of A. In [Appendix A we will prove that our 
results are independent on the specific choice of the maximally entangled input state of 
AB. Henceforth we will set the input state in the singlet state: 

Pab = |*-)ab(^-|, \'^-)ab = (|0)a® |1)b- |1)a® |0)b)/v^. (24) 



3.1. Collecting data in reflection and in transmission 

Let us focus first on the experimental setting in which the observer collects both 
transmitted and reflected signals of the probe A. 
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From the expression (17) of the QFI matrix it immediately follows that the reflection 
and transmission components of the state p*^{v) defined in Eq. (13) provide two 



separate contributions, as they are associated to orthogonal subspaces. Also, as in the 
case of the direct estimation discussed in the previous section, it turns out that the QFI 
matrix possesses a simple form in polar coordinates, where it is diagonal independently 



of the initial state of X. Indeed, upon diagonalization of the state (13) we find that in 



Cov[{r,0,<P)]>-H'^+^ir,e,<py 



this case the inequality (21) gets replaced by 

1 

M 
where the QFI matrix is 

with 

8^2(1 + 18^2 ^ g3^4 



t+r t+r t+r ■ 2 

r ,V 'V sm I 



.t+r. 



J+r, 



r,n) 

r,Q) 



(1 - r2)(l + fi2)(l + 5^2)^1 + 9^2)2 ^ 

1 



(25) 



(26) 



(27) 



(l + fi2)(l+9^2) 

X [4(1 + 5fi2)(l + 9fi2)(l + 18^2 _^ q3q4^) 

4^2 + 68fi^ + 720fi^ + 18631]^)] 
32r2fi2 



-r^il 



X 



[4(1 + 9^2)2 _ ^2(1 _ 9fi2)2][4(l + 5^2)2 _ ^2^^ + 3^2)2] 



f28) 



m 



Appendix B we also report the expression of the bound in cartesian coordinates). As 



in the case of the direct estimation, H^p^ does not depend upon the azimuthal angle 0, 
while it is a function of the radial coordinate r and the polar angle 9. Such a behavior is 
associated with the symmetry of the coupling Hamiltonian ([2]), which does not possess 
a preferred spatial direction, and of the input state of v4i?j^ It is also worth pointing 
out that the matrix H^ vanishes when Vt is zero or infinite (that is, infinite or zero 
incident momentum hk of the probe A). This implies that for such configurations no 
recovering of information on X is possible. Indeed if fc = it means that A is initially 
at rest and will not be able to be scattered by X, while for very large k it means that A 
propagates so fast along the line that the interaction with X can only have a minimal 
(asymptotically vanishing) impact on its evolution. 

As a specific example suppose then that the observer, already knowing the value 
of the parameters 9 and 0, is interested in recovering the missing parameter r, which 
determines the purity of the target system X, Ti p\{v) = r^. Equation (25) then yields 

1 1 1 _ ^2 



Var[r] > 



;i + fi")(l + 5fi^)(l + 9fi 



2\2 



(29) 



which should be compared with the quantity {l/M){l/cf^'^) = (1 — r'^)/M one obtains 



in the direct estimation case, i.e., Eqs. (21)-(23). Since the rational function of ^2 on 



^ We stress that the dependence of Hf~^ 



EA 



on 9 has nothing to do with the probe-target couphng: as 



in the case of H it is a geometric artifact of the polar representation. 
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(1 - r2)cf+^ 0.5 




Figure 3. Functional dependence of tlie rescaled QFI function (1 — r^)c*"'"'^, whose 



inverse bounds the accuracy ( 29 ) achievable in the determination of r with an EA 



strategy where AB is initialized in a maximally entangled state and data are collected 
both in transmission and in reflection cases. The dot indicates the optimal incoming 
momentum for probe A. 



the rhs of Eq. (29) is larger than 1, it follows that (l/M)(l/c^"") is always smaller than 



the EA bound (this is very much expected since in the case of direct estimation the 
observer has access to the X system, while in the EA strategy he can only recover 



info on X through the probe A). The minimum of the rhs Eq. (29) is reached when 
Q ~ 0.616 as shown in Fig. [3} this value provides the optimal input momentum hk of 
A via Eq. (|6|, for recovering the parameter r. For such a choice the EA strategy misses 
the direct estimation accuracy just by a factor of 1.52. Notice finally that for both 
the EA strategy and for the direct one, the accuracy bounds vanish with the a/1 — r^ 
distance from the surface of the Bloch sphere (this is a consequence of the fact that it 
is intrinsically simpler to distinguish pure states from mixed states). 

Consider next the case in which the observer, already knowing the value of the 
parameters r and 6, is interested in recovering the azimuthal phase (p of the target 



system X. This is given by the third diagonal element of the QFI matrix (26), i.e., 

1 



Var[0] > 



1 



(30) 



Mc*+"(r,fi)sin2^' 

which again should be compared with the bound (l/M)(l/cg 
(l/M)(l/r^ sin^ 6') one gets for the direct estimation. Again one can verify that the 



.dir g-j^2 I 



rhs of Eq. (30) is always larger than the direct estimation threshold (notice however 
that both expressions diverge for ^ = and 9 = tt due to fact that for such choices the 
variable is not even defined). For the sake of simplicity we focus on the special case 
in which X is in a pure state on the equatorial plane of the Bloch sphere (i.e., 6 = 7i/2 
and r = 1). In this case Eq. (30) yields 



Var[0] > 



1 3(1 + f^2)(l + 3fi2)(^ ^ 7^2)^^ ^ g^2) 



(31) 



M 32^2(1 + 10^2 + 27^4) 

which reaches its minimum when Q ~ 0.637, where the rhs is ^ 1.354/M, falling short 
of the direct threshold 1/M by 35%. 
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(1 - ■P)c* 




Figure 4. (Color online) Plots of the rescaled QFI coefficient (1 — r'^)c^{r,n) of 



Eq. (331 for the transmission case as a function of r and Q. For every value of r there 



exists an optimal value of the incident momentum k. The upper and lower curves in 
the right panel refer to the case of completely mixed (i.e., r = 0) and almost pure 
states (i.e., r ~ 1), respectively. 



3.2. Collecting only reflected or transmitted data 

By repeating the same analysis for the case in which only transmitted or reflected probes 



are detected (see Eqs. (14) and (15)) we find that Eq. (25) still holds with the QFI matrix 



H^j^{r, 9, 0) being replaced by 



H'Jj:{r,e, 



diag{cl^'',c*J'',cJ''sin^e), 



with 



(r,fi) 



{r,Q) 



2fi2[3(i + 3fi2)(^^ ^ 7g^2 ^ 117^4) _ 2^2(9 ^ 5q^2 ^ 45^4^] 
(1 - r2)(l + n2)(i + 5^2)(i + 9fi2)[9(i + 3^2)2 _ 4^2] ■ 



2\21 



3(1 + fi2)(l + 3fi2)(l + 9fi2)[4(l + 5f]2)2 _ ^2(1 + 3^2)2] ' 

in transmission, while 

2(]2[(i + 7fi2)(^ ^ 3g^2 ^ 207^^) - 2r^Q\l + 18^^ + 117^^^] 



c:{r,n) 



{r,Q) 



4r2fi2[2(i 



r2)(l + fi2)(l + 9fi2)2[(l + 7^2)2 _ 4^2^4] 

9fi2)(i + 36f]2 + 2071]4) - r2f]2(i _ 9^2)2] 



(32) 

(33) 
(34) 

(35) 
(36) 



(1 + fi2)(l + 7fi2)(l + 9fi2)[4(l + 9^2)2 _ ^2(1 _ 9^2)2] ' 

in reflection. 

Notice that, in this case, even if there is still a quadratic divergence of c^ and c* 
at r = 1, their dependence on r does not factorize. Hence in the estimation of the 
radius r, the bound (29) is replaced by a lower bound, (l/M)(l/cr ), which for every 
value of r admits an optimal momentum for the incident probe A that maximizes the 
associated QFI, as shown in Figs. |4] and |5} Interestingly enough, however, there exist 
"optimality intervals" for the incident momentum (i.e., Q G [0.51, 0.55] for transmission, 
and Q G [0.67,0.68] for reflection), which guarantee rather high performances for all 
values of r. Observe that in order to achieve the optimal estimation from transmitted 
data we have to send the probe faster than the case of reflection (recall that Q oc 1/k): 
this is in perfect agreement with the intuitive idea that to efficiently collect the data in 
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Figure 5. (Color online) As in Fig. l4J when only reflected data are collected. 



C, 3 




Figure 6. (Color online) Plot of the quantum CR bound for rescaled variance 
MVar[r], of the estimation of the radius r of the Bloch sphere describing the target 
qubit, in polar coordinates. The dashed line refers to the ultimate limit achievable 
by direct access to the target, while the solid lines refer to the case in which both 
reflected and transmitted, only transmitted and only reflected data are collected. In 
particular, the case of both reflected and transmitted data is plotted at the optimal 
value n ~ 0.61, which, in this case, is independent of the value of r. The case in 
which only reflected or only transmitted data are collected has been plotted instead by 
considering, for each given r, the corresponding optimal value of il. Notice that for all 
values of r, we find c*"'"'' > Cr , since by collecting all the scattering data we gain the 
largest amount of information on the target. Futhermore, observe that transmission 
and refiection give almost the same result, with a little improvement of transmission 
with respect to reflection. 



transmission A should be given a sufficiently large initial momentum to prevent back 
scattering from X (and vice versa for reflection). Between the above two intervals lies 
the optimal value of Q for the case in which we collect all scattering data, Q ~ 0.61. 
In Fig. [6] we plot the QFI of the radius r (purity) for the case of transmitted, reflected, 
or both reflected and transmitted data. Of course, the last case leads to the best 
estimation of the target purity, as we are collecting the largest amount of information 
from the scattering process. 
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Analogous considerations apply also to the estimation of the azimuthal angle (j). One 



can easily extract them by a direct comparison of the coefficients (34) and (36) with the 



corresponding expression associated with the transmission and reflection strategy given 



in Eq. (28) and also with the direct estimation procedure [see Eq. (23) 



4. Comparing EA and NEA strategies 

In this section we will present a comparison between the EA strategies introduced in 
the previous section with the NEA strategies obtained by restricting the analysis to the 
case in which the initial state pab of the probe A and the ancilla B is separable. For 
the sake of simplicity, we flnd it instructive to restrict the analysis to the situation in 
which, as in Fig. |6} the observer aims only to estimate the purity of the target state 
X. Speciflcally we will work under the assumption that the three dimensional vector v 
which specifles pxiv) in the Bloch sphere lies on the z axis, that is Vx = Vy = Q and 

PxM = ^^^, ^.€[-1,1]. (37) 



Under this condition the quantum CR bound (16) reduces to an inequality for the 
variance Var[t'2] of the unique parameter Vz, i.e., 

Var[..] > l^Hl\^/^il{vz)-\ (38) 

where HyJL-p'Avz) are the QFI functions of the problem computed as usual by 



'EA/NEA 



exploiting the spectral decomposition of the output states AB/A and using Eq. (17) 
(which in this case contains only derivatives with respect of the unique parameter Vz)- 
For the EA conflguration the functions H-^^ v^z) coincide with the third diagonal 
elements of the QFI matrices in cartesian coordinates evaluated at f ^ = fy = 0, namely, 

rrt+rf ._ 8^^(1 + 18^^ + 63^^) 

^^ ^^' (l-z;2)(l + fi2)(i + 5Q2)(i + 9^2)2' l-^y^ 



(1 - vl){l + Q?){1 + 5fi2)(i + 9fi2)[9(i + 3(]2)2 _ 4^ 
2^2[(i + 7^2)(^ ^ 3g^2 ^ 207^^) - 2vl^'^{l + 18^^ ^ 117^4^] 

''^^'""' ~ (1 - VDH + fi2)(l + 9^2)2 [(1 + 7^2)2 _ 4^2^4] ' ^ ) 



Notice that the above expressions can be retrieved from Eqs. (27), (33) and (35) 



by substituting r with Vz- The general expression of the QFI matrix for the EA 



configuration can be found in [Appendix B 



Before computing the QFI functions for the NEA configuration, it is worth 
observing that due to the convexity of the QFI function (e.g., see Ref. [16]), it follows 
that in the NEA scenario the contribution of the ancilla B in the scattering can be 
completely neglected, see Fig. |2] Indeed, since the maximum of a QFI function is 
always achieved on pure states, one can restrict the analysis of the NEA configuration 
to pure separable states of AB: for them however the output states of AB also factorize 
and one can completely ignore the subsystem B. Consequently, the states of the system 
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Figure 7. Envelope of the QFI for the NEA configuration corresponding to the optimal 
choice of 9a and fl for each value of v^ , if only transmitted data are collected. Notice 
the divergence at Vz = ±1, which for all il corresponds to 9a — and tt. 



after the scattering are given by Eqs. (|9]), (10) or (12), for the case in which we collect 
either all scattering data, or only transmitted/reflected probes. In particular we will 
assume 

1 + n- (Ta 



Pa 



n = (sin 9a, 0, cos 9a) 



(42) 



where, without loss of generality, the input azimuthal angle (pA has been set to zero by 
using the symmetry of the coupling Hamiltonian. 

Let us then consider first the accuracy achievable when collecting data only on 
transmission, i.e., assuming as output state the one given in Eq. (10). The resulting 
QFI is the following (involved) function of 9a'- 



HtqEAi^z,9A,^) 



[4(11 + 96^2 + ISm^) - vl{l + Q^){1 + 331]2) 

- 4:v,{8 + 43^2 + 3fi^) cos9a - 4(1 - fi^ + 8t;^fi^)(l + fi^) cos2^a 
- 4v,n\l + fi2) cos39a + v1{l + n'^y cos A9a] 



4(1 + 5f]2) _ t,2(i + 17^2) _ y^^i + n'^)(4cos9A - V, cos29a) 

fi2 



(1 + fi2)(3 + 9^2 _ 2v^ cos9a){1 + 7^^ + 2^,^^ ^Qg^^) 
Still one recognizes some general features we already observed in the previous section. 
In particular the QFI function vanishes at fi = and i7 = oo. Notice also that when the 
target state is pure, there exists a particular choice of the initial state of the probe such 
that the above quantity diverges quadratically, as found for the EA strategy, namely 
f^ = 1,6'^ = and Vz = —1,9a = ti", and arbitrary (finite) i7 > 0. Finally, notice that 
for each value of Vz there exists an optimal point for the pair {9a, ^) that maximizes the 
QFI. In Fig. It] we plot the envelope of the QFI for Vz G [—1, 1], that is the maximum 
value of H^^j^{vz), which can be proved to be a symmetric function of Vz- Notice that 
the initial state of A that yields the best estimation of the target state depends on 
Vz, that is, on the initial state of the target, which is in principle unknown. We can 
only say that once we have set the direction of the probe A in the Bloch sphere before 



(43) 
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Figure 8. (Color online) Envelope of the QFI for the NEA strategy corresponding to 
the optimal choice of Oj^ and il for each value of v^, for the case in which we collect 
either transmitted, or reflected or all scattered probes. Notice that the QFI diverges 
for Vz — l^ and —1 when 9a = and tt, respectively. 



the scattering, the best optimization we can get involves pure target states with Bloch 
vectors parallel to this direction. The same analysis can be repeated for reflection and 
yields: 

= [4(5 + 23^2) _ t;2(l + ^2^) _ 4^^(^3 _ 2^2^) ^^^q^ 

+ 4(1 - 5f]^) cos 29a - 4:V,{1 + 2n'^) cos 39a + vl{l + f]^) cos 4:9a] 

1 
^ 3(1 + 3fi2)(i + 7^2) _ 2t;2fi2 _ 2t,^(i + 4f]2 _ 9^4) ^os 9a - 2v^,n^ cos 29a 

X . (44) 

(l + fi2)(4_^2_4^^cos^A + ^^^cos2^^) ^ ^ 

Similarly to the case of transmission the optimal values of 9a and Q are symmetric 
functions of Vz- Moreover, when the target state is pure, i.e., Vz = 1 and —1, the QFI 
diverges for 9a = ^ and tt, respectively, as for transmission. 

Finally, the largest amount of information on the system can be inferred by 
collecting both transmitted and reflected data. In this case the QFI associated to the 
state ([9]) of the incident probe after the scattering is given by 

H'^^l{vz.9A,n) 

= {{l + 7n^ + 2vzn^cos9A) 

X [4(1 + 5^2) - vl{l + 17^2) - Vzil + n^)i4:Cos9A - VzCos29a)] 

X [4(5 + 27fi^) -vl + A{l- 9f]2) cos 29a - I6vz cos 9a^ + vl cos A9a] 

+ (4 - i;2 - Avz cos 9a + vl cos 29a) [3(1 + 3^'^) - 2vz cos 9a] 

X [4(3 + 48^2 + 181^4) - vlVl?{l + 33^2) - l2vzVl?{l + fi^) cos 9a 

- 4(1 + 8fi2 - fi^ + Svln^) cos29a - w^fi^(l + fi^)(4cos3^A - VzCos4:9a)]} 

1 
^ 4(1 + 5fi2) - t,2(i + 17^2) _ ^^(1 + n'^){Acos9A - Vz cos29a) 
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1 



(4 - i;2 - Av^ cos Oa + vl cos 26 a) (3 + 9^2 _ 2v^ cos 6 a) 

This function exhibits the same properties as H^^p^ and H^^a- The three curves plotted 
in Fig. [8} all symmetric with respect to Vz, refer to the optimal values of the QFI for the 
three cases analyzed above. As expected, the collection of all scattering data returns 
the best tomographic reconstruction of the target state. Notice that the transmission 
strategy seems to overcome the reflection one, like in the EA strategy. Also it can be 
shown that for all Vz the optimal value of Q for the transmission is lower than that for 
reflection, and for the case in which we collect both transmitted and reflected data it is 
in between (this results have been obtained by numerical optimizations). 

In order to compare the efficiency of the NEA strategies with the EA strategies 
discussed in the previous section, in Fig. [9] we plot the maximum QFI for fz > 0, for the 
EA (solid line) and the NEA (dashed line) strategies (the expression for the EA strategies 
have been obtained by exploiting the representation of the QFI matrix in cartesian 



coordinates reported in Appendix B). The role played by the entanglement between the 
probe and the ancilla before the scattering is evident: it implies an enhancement in the 
QFI for all Vz, both for the case in which we collect all the scattering data and for the 
case in which we have access only to trasmission/refiection events. 



5. Conclusion 

We have presented a detailed study of the tomographic state reconstruction of a target 
system, obtained by monitoring a scattered probe. Focusing on the special case in which 
both the target and the probe are qubit systems, and assuming the scattering to take 
place on a 1-D line, we used quantum estimation techniques to evaluate the efficiency 
of the process in several configurations of interest. In particular we distinguished two 
regimes: the EA regime in which the observer is allowed to initialize the probe in an 
entangled state with an external ancilla that it is kept in the laboratory; and the NEA 
regime where instead no entanglement is allowed between the probe and the ancilla. As 
expected, when all the other settings are kept identical, the EA strategies turn out to 
be more effective then their NEA counterparts. 

Within both regimes we have also studied what happens when the observers have 
access to all or only part of the scattered data. Specifically we consider the cases in which 
only transmitted or refiected data are used in the tomographic reconstruction, noticing 
that these regimes are characterized by different optimal values for the input momentum 
of the probing particle (the transmitted scenario being characterized by higher optimal 
input momenta than the refiected one). We have also analyzed the situation in which 
both transmitted and refiected data are available to the observer, assuming though that 
no joint coherent measurements could be performed on the associated quantum degree 
of freedom (i.e. we explicitly excluded the possibility of performing joint detection on 
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Figure 9. (Color online) Maximum value of the QFI associated to the parameter Vz 
of the target qubit X. The solid line refers to the EA strategy, where the probe A is 
in a maximally entangled state with the ancilla B, while the dashed line refers to the 

NEA. 



the left and right side of the 1-D channel, an hypothesis which is very much reasonable 
if the detectors are located at sufficient large distances from each other). The overall 
accuracy clearly benefits from this possibility: still it remains below the threshold [3] 
one gets when direct access to the target system is allowed. An open problem is to the 
determine whether or not one could exploit other sort of quantum resources to close 
such gap. A possible candidate which we aim to explore in a future development of 
the work, is to allow the observer to entangle different probes together (in the present 
scenario indeed, even though we used M probes, they were all prepared in a factorized 
configuration of the same initial state). This strategy could in fact benefit from super- 
additivity effects arising from the presence of non classical correlations between the 
various probes, resulting in higher performances of the tomographic reconstruction. An 
interesting open problem is also to understand to what extend the specific features we 
have observed in our simple scattering model (qubits interacting along a 1-D line via 
Heisenberg-like coupling) could be generalize to more complex configurations (say by 
increasing the size the of the target and/or of the probe, or by allowing the scattering 
to take place in a 2-D or a 3-D setting). 

This work was supported by the Italian Ministry of University and Research through 
FIRB-IDEAS Project No. RBID08B3FM as well as under the bilateral Italian- Japanese 
Projects II04C1AF4E on "Quantum Information, Computation and Communication," 



Entanglement- assisted tomography of a quantum target 19 

by the Joint Italian- Japanese Laboratory on "Quantum Information and Computation" 
of the Italian Ministry for Foreign Affairs, and by a Special Coordination Fund 
for Promoting Science and Technology and a Grant-in-Aid for Young Scientists (B) 
(No. 21740294) both from the Ministry of Education, Culture, Sports, Science and 
Technology, Japan. 

Appendix A. Equivalence between maximally entangled states 

We explicitly prove that our analysis for the EA strategies is independent of the 
particular choice of a maximally entangled state for the subsystem AB. 

The Heisenb erg-type coupling ax ■ <ta describing the interaction between qubits X 
and A can be written in terms of the so-called swap operator Sx\a as 

(Tx-(Ta = 2Sx\a-^xa, <Sx\A = '^\i)x{j\®\j)A{i\- (A.l) 

It is a unitary self-adjoint operator and is characterized by the following very simple 
property 

(Ix ® Ua)Sx\a(Jx ® U\) = (f/i ® Ia)S'x\a{Ux ® W, (A.2) 

where Ux and Ua act in the same way on Tix and T-La, respectively {Ux = Ua), and 
5^,^ is the swap operator in the rotated frame 

<Sx\A = X^ \ri)x{rj\ ® \rj)A{ri\, \ri)x/A = Ux/A\i)x/A- (A. 3) 

The above property is trivially conserved for the interaction Hamiltonian crx ■ cta- The 
equivalence among all the maximally entangled input states for the EA strategies is 
thus straightforward. Indeed, a generic maximally entangled state of the probe A and 
the ancilla B can always be written as \^uv)ab = {Ua ® VbiY^') ab-, with \^~)ab the 



singlet state (24). If we send \^uv) ab to the target qubit X it can be easily shown that 
the contributions given by Trx{. . .} in the final state of the subsystem AB become 

{U\ ® Vl)i:ix[{UAStAUA^){px{v) ® \m-)AB{'^-mASfAU\)]{UA ® Vb) 

= {U\ ® V},)1:ix[S'"'xa{Px{v) ® \^-)AB{'^-\)SfA\UA ® Vb) 
with 

SYa = «t,r(fi) + /3t,r(^)(25'x|A - Ixa), p'x = UxPxUJ,. 

Since the physical properties of the system do not depend on the choice of the reference 
basis, the results of our analysis for the EA strategies are completely independent of the 
choice of an initial maximally entangled state of AB. 

Appendix B. QFI matrix in cartesian coordinates 

In this appendix we consider the explicit expression of the Fisher matrix for the EA 
strategy and discuss the symmetry properties of its elements. 
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For the case in which both transmitted and reflected probes are collected, the off 
diagonal elements of the Fisher matrix are given by 

[-"EA {'^x,Vy,Vz)\ij,iy^j 



2viViil'^ 



(1 - V2)(l + fi2)(l + 5^2)(1 ^ 9^2)2 

X {(1 + 7fi2)(l + 9fi')(3 + 131]')'[4(1 + 9ny - v\l - 9n^f] 

+ 3(1 + 3fi2)(i + 5^2^^^ ^ 271]')2[4(1 + 5nY - «'(1 + 3^1^]} 

1 

^ [4(1 + 9f]2)2 _ ^2(1 _ 9^2)2] [4(1 + 5^2)2 _ t;2(l + 3^2)2] ' 

with i,j G {x,!/, z}, while for the diagonal elements we have 

I. EA {"^xyVyyVzJlii 

2D? 



(B.i; 



(1 - V2)(l + (]2)(1 + 5^2)(1 + 9^2)2 

X {(1 + 9f]2)(3 + 13fi2)[4(l + 9^2)2 _ ^2(1 _ 9^2)2] 

X [t;,'(l + 7f]2)(3 + 13fi') + 4(1 - v^){l + 5^')=^] 
+ (1 + 51]2)(1 + 27f]2)[4(l + 5^2)2 _ ^2(1 ^ 3^2^2j 

X [?,vl{l + 3fi2)(l + 27f]2) + 4(1 - v2)(i + 9f]2)2j| 

^ [4(1 + 9fi2)2 _ ^2(1 _ 9^2)2] [4(1 + 5^2)2 _ ^2(1 + 3^2)2] ' (^-2) 

with i,j G {x,y,z}. Notice that if there is only one out of the three parameters 
characterizing the initial state of the target different from zero, f « 7^ and Vj = f ^ = 0, 
the QFI matrix is diagonal and thus its ii element coincides with the single parameter 
QFI for f j. Furthermore, due to the symmetry of the Heinsenberg-type interaction and 
of the singlet state of A and B before the scattering, we find 

(l-v^)[Hl\^{v,,Vy,v,)],, = (1 -r2)4+^ I G {x,y,z}. (B.3) 

If we are able to detect only transmitted or reflected probes we get 

[-"EAl'^a;) "^y; ''^z)\ij,i^j 

2ViVjV? 



3(1 - t;2)(i + f]2)(i + 3^2)(i + 5^2)(i + 9^2) 

X {[3(1 + 31]2)(i + 7l]2)(3 + i3f]2)2[9(i ^ 3^2)2 _ 4^2^ 
+ 8(1 - v^){l + 5f]')[4(l + 5fi2)2 _ ^2(1 ^ 3^2)2]| 

^ [9(1 + 3^2)2 _ 4^2] [4(1 + 5^2)2 _ ^2(1 + 3^2)2] ' ^^'^^ 



[-"EAl'^^a;; Vy, Vz)\ii 

2^2 



3(1 - V2)(l + fi2)(l + 3(^2)(1 + 5^2)(1 + 9^2) 

X {2(1 - V^){1 + 5fi2)[4(l + 5^2)2 _ ^2(1 ^ 3^2)2| 
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X [9(1 + 3nY - A{v] + vl)\ 
+ 3(1 + 3fi^)(3 + 13(]2)[9(1 + "iVL^f - 4i;2] 

1 

[9(1 + 3fi2)2 _ 4^2] [4(1"^ 5^2)2 _ ^2(1 + 3^2) 

and 



^ rn/1 , or.2\2 /f„.2ir^/i , c:r.2\2 ..2^ , or.2N21' (^•^) 



[-"EAl'^a::' Vyi Vz)\ij,ijLj 



2ViVjVL^ 



(1 - t;2)(l + fi2)(l + 7^2)(l + 9^2)2 

X {3(1 + 3fi2)(i + m''){l + 27fi2)2[(i + 7^2)2 _ 

- 8(1 - v^)Q!^{1 + 9n2)[4(l + 9^2)2 - ^2(1 

1 


- 9Qy 


^^ (1 + 7^2)2 _ 4i;2^4][4(l + 9^2)2 _ ^2(1 _ 9^ 

J^EAK'^XJ Vy, Vz)\ii 

2^2 


2-\21' 



(B.6) 



(1 - V2)(l + fi2)(l + 7^2)(1 + 9^2)2 

X {(1 + 71]2)(1 + 271]2)[(1 + 7^2)2 _ 4^2^4] 

+ 2(1 - t;2)fi2(l + 9f]2)[4(i ^ 9^2)2 _ ^2(^ _ 9^2)2| 

x((l + 7fiY-4fi'(^J-^D)} 

^ [(1 + 7^2)2 _ 4i;2^4][4(l + 9^2)2 _ ^2(1 _ 9^2)2] ' (^•'^) 

Vi, j, A; G {x, y, 2;}. Analogously to the case in which we detect all the scattered probes, 
we find that also for the transmission and the reflection cases if we have f , 7^ and 
Vj = Vk = 0, i y^ j j^ k E {x,y,z}, the QFI for t>j is diagonal, and furthermore an 



equation similar to (B.3) holds: 



(1 - V,') [ifEA (^x, Vy, v^)U = (1 - r2)4/^ z G {x, y, z}. (B.8) 
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